Quantum dynamics and breakdown of classical realism in nonlinear oscillators 
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The dynamics of a quantum nonlinear oscillator is studied in terms of its quasi-flow, a dynamical 
mapping of the classical phase plane that represents the time-evolution of the quantum observables. 
Explicit expressions are derived for the deformation of the classical flow by the quantum nonlinearity 
in the semiclassical limit. The breakdown of the classical trajectories under the quantum nonlinear 
dynamics is quantified by the mismatch of the quasi-flow carried by different observables. It is shown 
that the failure of classical realism can give rise to a dynamical violation of Bell's inequalities. 
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Introduction Single-degree of freedom nonlinear 
quantum oscillators are familiar as simple model sys- 
tems and textbook examples of quantum mechanics. 
They are also of direct physical interest as effective 
Hamiltonians for the pardigm systems of trapped atoms 
and ions [3, 0, S] and single-mode cavity quantum 
electrodynamics and quantum optics 0, d, Q . 

These systems are often studied experimentally by ex- 
citing them with an external perturbation and following 
the resultant evolution. The excited states are in general 
completely arbitrary: They are not necessarily close to 
specific energy states, nor to a superposition of a small 
number of energy states. For example, a nonlinear single- 
mode cavity can be excited into an arbitrary superposi- 
tion of photon states, and a trapped atom may be excited 
into a state with an arbitrary wave function and in 
principle these states can be tailored for specific purposes 
[7( . At this stage it becomes important to develop a fun- 
damental understanding and practical calculational tools 
for quantum nonlinear dynamics of oscillators of arbi- 
trary form, with an arbitrary initial state. In this way, 
for example, it would be possible to predict the behavior 
of a nonlinear cavity as a light source. 

In this Letter we put forward a method of general ap- 
plicability for the study of the quantum dynamics of non- 
linear oscillators, derived from the classical nonlinear dy- 
namics. Since the classical dynamics is simple and well- 
understood, it serves as an excellent starting point for the 
analysis of the quantum dynamics. In this manner, the 
time evolution of quantum observables is studied here for 
oscillators of arbitrary form and strength of nonlinearity, 
with an arbitrary initial wavefunction. It is obtained as 
a power series in Planck's constant, where the semiclas- 
sical limit is approached when the level spacing changes 
slightly between one pair of levels and the next one, an 
assumption which is certainly justified in the applications 
cited above. On the other hand, it is not assumed that 
the uncertainties are small, nor is any other assumption 
made on the properties of the quantum state. 

The dynamics of the oscillator is studied via its quasi- 
flow, a time-dependent mapping of the classical phase 
plane to itself that yields the dynamics of the observ- 



able for a given initial wavefunction. The quasi-flow is a 
smooth function of h which reduces to the classical flow 
when h — > 0. The analysis of the dynamics via the quasi- 
flow offers several advantages as compared with standard 
semiclassical methods. Unlike the linearization method 
[To| no assumptions are made on the smallness of the 
uncertainties. We do not rely on the calculation of highly 
oscillatory i nteg rals needed in the theory of semiclassical 
propagator [ill Il2| , and in this manner avoid caustic 
singularities. Furthermore, the analysis is exact in the 
nonlinearity, facilitating the explicit calculation of the 
quasi-flow for a large class of oscillators without recourse 
to numerical simulations, enabling the quantitative study 
of important fundamental effects arising from the quan- 
tum nonlinearity. 

Several nonclassical effects are identified and analyzed 
in terms of the quasi-flow. Squeezing [f| Q in its gen- 
eral sense refers to the distortion of the uncertainties by 
the classical nonlinearity, which sometimes produces sub- 
ground state fluctuations in one of the quadratures; in the 
present analysis squeezing is obtained almost trivially, 
since it occurs also for classical phase-space distribution 
undergoing nonlinear dynamics. A second effect is the 
deformation of the quasi-flow by the quantum nonlinear 
dynamics. In the semiclassical domain the quasi-flow de- 
formation manifests as a drift in the mean values of ob- 
servables. In particular, energy conservation combined 
with the noncommutativity of the observables implies a 
quantum drift in the amplitude of the oscillations. 

Of course, quantum mechanics cannot be fully de- 
scribed in terms of a phase-space flow. This basic fact 
is reflected in the present analysis in that the quasi-flow 
is observable-specific. That is, each observable carries a 
quasi-flow modification, and the mismatch of quasi-flows 
is responsible for a fuzziness which spoils the classical 
realism that may be attributed to the initial state. We 
demonstrate this dynamical breakdown of classical real- 
ism by calculating explicitly the quantum discrepancy for 
a generating function of observables, and showing that it 
can lead to a dynamical violation of Bell's inequalities. 

Classical dynamics of integrable systems is most conve- 
niently traced using action-angle variables. The dynam- 
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ics in the original variables is then expressed as a compo- 
sition of three canonical transformations. The quantum 
analog of the action-angle transformation is the unitary 
transformation to quantum normal form, developed in 
Refs. ID, 14 1 . Accordingly, we first calculate the quasi- 
flow of oscillators in normal form, which describes the 
dynamics of nonlinear oscillators up to unitary equiv- 
alence. Then we derive a composition formula for the 
phase-space representations of unitary transformations, 
thereby implicitly obtaining the quasi-flow for nonlinear 
oscillators in general form. 

The nonlinear oscillators studied here are ideal system, 
where interactions with the environment and decoher- 
ence are neglected. These important effects lend them- 
selves naturally to the dynamical approach taken here, 
and their study is planned in future works 

The Weyl quasi-flow As is well known, although po- 
sition and momentum are not simultaneously measurable 
in quantum mechanics, it is still possible to define phase- 
space quasi-probabilites [17[. The quasi-probabilities 



is, there exists a unitary operator U such that 



have to be accompanied by a symbol map, which assigns 
a phase-space function to each Hilbert space operator, 
such that the expectation value of the operator is equal 
to the phase space integral of the symbol with the quasi- 
probability as a measure. 

The quasi-flow is the time-dependent mapping z t (z) of 
phase space to itself defined by the symbols {q t {z) , Pt{z)) 
of the canonical position q t and momentum p t in the 
Heisenberg picture (we use the symplectic notation z 1 = 
—Z2 = q, z 2 = z\ = p with implied summation over re- 
peated indices). Since the symbols of q and p are q(z) = q 
and p(z) = p (respectively) the initial value of the quasi- 
flow is the identity mapping. The quasi-flow facilitates 
the calculation of expectation values of the canonical 
variables as a function of time for any initial state, and of 
any other observable by means of a static transformation, 
as demonstrated below. 

Evidently, the form of the quasi-flow depends on the 
choice of the quasi-probability and symbol scheme. We 
choose here the Wigner-Weyl scheme [T3], since it has 
several convenient properties, among them the simple 
form of the semiclassical expansion of the symbol of the 
product (the Weyl-Groenewold-Moyal product rule [l^|). 
The quasi-flow may be defined for any value of ft, but it 
has a particular advantage of reducing to the classical 
flow in the limit h — > 0, enabling us to take full advan- 
tage of the simplicity of the classical flow. 

The quasi-flow of nonlinear oscillators in normal form 
Nonlinear oscillators are defined here as (spinless) one- 
degree of freedom quantum dynamical systems governed 
by Hamiltonians whose Weyl symbol, the classical Hamil- 
tonian H(q,p), has a single nondegenerate phase space 
minimum and no other fixed points. It has been demon- 
strated in Ref. [HI that such Hamiltonians are unitar- 
ily equivalent to a normal form Hamiltonian which com- 
mutes with the harmonic oscillator Hamiltonian. That 



Q = uHu, 

H = f{B) , 



P = U^pU , 

z_ Q 2 + P 2 
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The classical limit of this transformation is a Cartesian 
version of the canonical transformation to action-angle 
variables; in particular the classical limit of B, the symbol 
of B, is the classical action. Like its classical analog, the 
quantum normal form Hamiltonian is a representative of 
its equivalence class where dynamics takes a particularly 
simple form. For this reason we will assume here that the 
Hamiltonian has been reduced to normal form, construct 
its quasi-flow and study its properties, returning to the 
general case in the last part of the paper. 

It is convenient to formulate the normal form dynamics 
in terms of the complex coordinate f3 = -^(Q + iP) an d 

the annihilation operator b — ;4j(Q + iP) (note that 

b differs by a factor of \fh from the usual definition). 
The symbol of bt will be denoted by /?*. In complex 
coordinates B = |/3| 2 and B = \{btf + tfb). 

Using the algebraic properties of the observables [lij ]. 
their dynamics in normal form can be expressed as a 
rotation with an operator-valued frequency 



(2) 



where U±{B) = ±\{f{B ± h) - f{B)). Note that 
lim Q,± = u) is the classical (action-dependent) frequency. 

Either of the two expressions on the right-hand-side 
of the last equation forms a useful starting point for a 
semiclassical expansion, since their factors have symbols 
with smooth dependence on H. This fact makes the calcu- 
lation of the semiclassical expansion of the normal-form 
quasi-flow a straightforward application of the rules of 
the Moyal algebra; the resulting expression to second or- 
der in h is 



P t (0)=ti c \l + tfft" +0{ff)) , 
where (3^ — fie^ 11 ^^ 1 is the classical flow, and 
Re/3 t (2) = \({u'{B)tf + BJ{B)J'{B)t 2 ) , 



(3) 



(4) 



Im/J 



(2) 



l f (-2B(w , (B)t) 3 + 5u"(B)t + 2Bu'"{B)i) , 

(5) 
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are the leading quantum correction to the modulus and 
phase of 0t, respectively. 

The physical significance of Eqs. (J3] E]) is as follows: 
The leading term describes the squeezing of the moments 
of b by the phase-plane shearing implied by the action- 
dependence of the classical frequency. Since the leading 
term in the dynamics of all observables is given by the 
classical flow, this effect is equivalent to the drift and 
deformation of the initial (Wigner) quasi-distribution by 
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the classical flow. The higher order terms, in contrast, 
express a nonlinear quantum drift, which has no classical 
analog. The imaginary part of the quantum drift gen- 
erates an accelerating drift in the phase of (b t ), and the 
real part generates a drift in the amplitude of the oscil- 
lations; at first glance the existence of an amplitude drift 
seems surprising, since energy conservation implies con- 
servation of the action operator B. In fact, it is precisely 
the conservation of the action, which generates the am- 
plitude drift through the Moyal algebra. Thus, on one 
hand B t (f3) — B, and on the other hand 



B + h 2 (2Bpi 2) - ^d lt dv$ B} &*ff'$ :) ) + 0(h 4 ) , (6) 



(c); 



where * stand for the Moyal product, <9 M = J^-, and 



The two results imply that Re (3\ 



(2) 



T ^ E d^d u f3{ c) d ll d v p ( t c) , consistently with Eq. flU). 

Dynamical breakdown of classical realism Unlike the 
leading (classical) term of the quasi- flow, the higher order 
terms cannot be directly used to derive the dynamics 
of the quasi-distribution. That is, it is not possible to 
calculate the Wigner function of a state at time t, given 
the initial Wigner function and the quasi-flow Eqs. ([3] [S} • 

Thus, in particular, the quasi-flow in itself does not suf- 
fice for the calculation of higher moments of the canonical 
observables, such as the uncertainties. For this purpose 
the quasi-flow has to be supplemented by the discrepancy 
dt{A) of an observable A at time t, being the difference 
between the symbol of A t and the symbol of A evalu- 
ated at the quasi-flow /3 t . The quantum discrepancy of a 
general observable in the q, p algebra is extractable from 
the discrepancy of the generating function e abt ~ aht of the 
complex parameter a, 



d t (e ab <- Qbt ) =tPe a % e) -*f>i B) 

x {Pxu'{B)t + P 2 {u'(B)t) 2 + P 3 (cj'(B)t) 3 

+ P A u"{B)t + P b uj\B)uj" {B)t 2 ) + 0(h 4 ) , (7) 

where P\ , . . . , P§ are explicit low-degree polynomials in 
B, \a\ 2 , aftl and a/3 t . Together with the quasi-flow, 
Eq. |(7J) gives the full information on the quantum dynam- 
ics of the (normal-form) nonlinear oscillator to 0(h 4 ). 

The irreducibility of the dynamics to a phase-space 
flow, as expressed by the quantum discrepancy, is a con- 
sequence of the breakdown of classical realism under the 
quantum nonlinearity. That is, suppose that the initial 
state is such that q and p have a classical hidden variable 
underpinning; the underpinning could be provided, for 
example, by the initial Wigner function, if it is nonneg- 
ative on the entire phase plane. Then the hidden vari- 
able underpinning breaks down once the system starts to 



evolve, and in general fails for q t and pt for any positive 
t, unless the dynamics is linear. 

It was pointed out in Ref. [l6| that the breakdown 
of the hidden variable underpinning can be detected 
when the classical symbol of the observable attains values 
which do not belong the spectrum of the observable as a 
Hilbert space operator. Following Ref. such observ- 
ables will be termed improper. When the hidden variable 
underpinning is provided by a nonnegative Wigner func- 
tion, its classical symbol is the Weyl symbol. It is a 
property of the Weyl calculus that the symbol of a func- 
tion f(S) of a linear combination S of Q and P is equal 
to f(S), and thus such observables are proper. However, 
if / is nonlinear, then the quantum discrepancy d t f(S) 
is nonzero for t > 0, and /(St) may become improper. 
That is, the quantum discrepancy is a direct measure of 
the breakdown of classical realism. As an example we 
may consider the observable C = cos aP (with symbol 
C), whose spectrum is bounded between —1 and 1. It 
follows from Eq. ([7]) that (in real coordinates) 

Ct(Q,o) = i+ft 2 ^y(?)+QV(^y(f)i 2 +o(fi 4 ) 

(8) 

which implies that for a large class of nonlinear oscil- 
lators, C t is improper for all t > 0, and the failure of 
classical realism increases quadratically in time. 

An important further consequence of the breakdown 
of classical realism arises when the nonlinear oscillator is 
entangled with other degrees of freedom. The entangle- 
ment is assumed to result from the initial preparation, 
after which the nonlinear oscillator evolves purely under 
the single-degree of freedom Hamiltonian. If the entan- 
gled state, like the EPR state 15|, 1(| , has a nonnegative 
Wigner function, then initially observations of Bell oper- 
ators constructed from functions of linear combination of 
Q and P are guaranteed to satisfy Bell's inequalities. At 
later times, however, measurements of a Bell operator 
constructed, for example, from cos(aQt) and cos(aPt) 
are liable to lead to violate Bell's inequalities, since these 
operators are improper for t > 0. As demonstrated by 
Eq. (JSj> , this nonlinear dynamical Bell inequality viola- 
tion increases quadratically in time. Further study of 
this topic is beyond the scope of the present work. 

Dynamics of nonlinear oscillators in general form. 
Since every nonlinear oscillator is unitarily equivalent to 
a normal form oscillator [l3| , the results presented above 
yield the quasi-flow of any given nonlinear oscillators af- 
ter it has been transformed to its normal form. Never- 
theless it is of interest to express the result in terms of 
the original observables q,p rather than the transformed 
observables Q, P. 

In the classical limit, the nonlinear oscillator flow can 
be expressed as a composition of three canonical transfor- 
mations (q,p) — ► (A, 4>) — ► (A, <f) + Lo(A)t) — ► (qt,Pt), 
where the first and last arrows are the transformation to 
action-angle variables and its inverse (respectively), and 
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the middle arrow is the time-evolution canonical trans- 
formation. Analogously the quantum evolution is the 
composition of three unitary transformations, (q,p) — ► 
(Q,P) — ► (Qt,Pt) — ► {Qt,Pt)- Each unitary trans- 
formation can be represented by its Weyl symbol. The 
Weyl representation for the normal form transformation 
was derived in Ref. (l3j , while the representation for the 
time evolution is the quasi- flow given by Eqs. ((3H5]). 

Unlike classical canonical transformations, however, 
the Weyl representation of the composition of the trans- 
formations is not equal to the composition of Weyl rep- 
resentations. Therefore, in order to obtain the general- 
form quasi-flow, a composition formula for Weyl symbols 
is needed. That is, given that the Weyl symbol of A 
with respect to the canonical set Q,P is A(z), and that 
the symbols of Q and P with respect to the set q,p are 



J 



Z^(z), to find the Weyl symbol a(z) of A with respect 
to q, p. In the following we derive a semiclassical expan- 
sion of a(z) to 0(h 4 ), which, together with the previous 
results yields the quasi-flow z^(z) of nonlinear oscillators 
in general form. 

The composition formula is derived directly from the 
definition of the Weyl symbol. By definition a(z) = 
TrAA(z), where A(z) = J* ^.e*(*" _ *"K are the basis 
operators for the Weyl representation [13] with respect 
to the set z, while A = J ^A(z)A(z), where A(z) are 
the basis operators with respect to Z. Therefore 



a{z) 



dz' 



A{z')TiA(z')A(z) 



J dz'A(z')5(z'] 



z) 



(9) 



where 8(z';z), the symbol of A(z') with respect to z divided by h, can be calculated with the help of the Moyal 
calculus, giving 

5{z'- z )= fi - » a ((^ w M *M - (^rt^z^nH^))^ - ^W) + ois 1 ) • (io) 

Using this result back in Eq. © gives the composition rule 
a(z) = A(z,Z(z)) , A(z,z') =(l- % 2 ({^d k d l Z^d k d l Z^d^ + {^d k Z x d k d l Z»diZ»)d x d^d v )\ A(z') . (11) 



Conclusions. The phase-space analysis of the dynam- 
ics in the Heisenberg picture was shown to furnish a 
powerful tool for the study of a broad class of nonlinear 
quantum systems. Fundamental nonclassical dynamical 
effects, including the dynamical breakdown of classical 
realism were studied nonperturbatively. The main re- 
sults are the explicit expressions for the quasi-flow and 
the quantum discrepancy in terms of the nonlinear fre- 
quency function to(B), which hold all the dynamical in- 
formation on the nonlinear oscillator and imply the non- 
classical phenomena. 
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